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We report measurements of dm/dn in Si MOSFET, where m is the magnetization of the two-dimensional 
electron gas and n is its density. We extended the density range of measurements from well in the metallic 
to deep in the insulating region. The paper discusses in detail the conditions under which this extension is 
justified, as well as the corrections one should make to extract dm/dn properly. At low temperatures, dm/dn 
was found to be strongly nonlinear already in weak magnetic fields, on a scale much smaller than the char- 
acteristic scales, expected for interacting two-dimensional electron gas. Surprisingly, this nonlinear behavior 
exists both in the dielectric, and in the metallic region. These observations, we believe, provide evidence for 
strong coupling of the itinerant and localized electrons in Si-MOSFET. 
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Magnetic properties of strongly-interacting electron 
|gas have long been a subject of intensive theoretical 
and experimental investigations. Coulomb interaction 
through the exchange term favors parallel spins, and 
therefore leads to ferromagnetism, whenever is strong 
enough. The strength of the interaction is determined 
by the ratio between the typical Coulomb and kinetic 
energies, and customarily characterized by the ratio be- 
tween the typical inter-electron distance p, p — l/s/nn 
in two dimensions, and the effective Bohr radius: r s = 
p/a*B IB- When the carrier density n decreases, the rela- 
tive effect of the interactions becomes stronger, contrary 
to the naive expectations. At sufficiently low density a 
clean system of itinerant electrons was predicted to be- 
come ferromagnetic, the phenomenon called Stoner in- 
stability [5] ; though in two dimensions it can occur only 
at zero temperature. 

A real system of itinerant electrons is always dis- 
ordered to some extent. Interactions, in the limiting 
case of strong disorder, favor antiparallel arrangements 
of the neighboring spins. The ground state of a pair 
of localized spins is a singlet, much like a Hydrogen 
molecule. A system of many localized spins preserves 
the tendency to antiferromagnetic order, see [3J for re- 
view. Intensive investigation of the magnetic properties 
of the doped semiconductors, particularly phosphorus 
doped Si in 80s [3J, lead to substantial understanding 
of the interplay between interactions and disorder. The 
observed divergency of the susceptibility with decreas- 
ing temperature at the low-density side of the metal- 
insulator transition was well understood; what remained 
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unanswered, is the divergent susceptibility at the high- 
density side of the transition. 

Interactions in 2D systems are more important than 
in 3D, which explains the wealth of experimental ob- 
servations and difficulty in their explanation. Unfor- 
tunately, extremely small number of electrons hinders 
thermodynamic measurements of 2D systems. On the 
other hand, two-dimensional gated structures, particu- 
larly Si-MOSFETs, provide unique possibility of gradu- 
ally changing the electron density, and, thus, interaction 
strength. Interest in the magnetic properties of 2DEG 
was sparked by the observation of strong suppression 
of conductivity in Si-MOSFETs by in-plane magnetic 
field [1]. The strong magnetoresistance was interpreted 
in [3 [6] as a quantum phase transition into a ferromag- 
netic state at the density of the metal-insulator transi- 
tion, n c . This interpretation was contested in [7J H] on 
the basis of Shubnikov-de Haas measurements, and in [9] 
on the basis of thermodynamic magnetization measure- 
ments. 

In the thermodynamic method, developed in [5] , 
recharging current, generated between the 2DEG and 
the gate in response to modulation of external mag- 
netic field B, is used to determine the derivative of 
the 2DEG chemical potential p2D with respect to the 
field. This derivative, by virtue of the Maxwell relation 
du/dB — —dm/dn, can be translated into the deriva- 
tive of the magnetic moment m at a given field with re- 
spect to the density n. One can, in principle, integrate 
—dm/dn over n to get m(n), assuming m is known at a 
certain density. In Refs. O HH] the integration constant 
was taken from transport measurements at some high 
density value. 
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It would be desirable to abandon this assumption 
by straightforward integration from n = 0, since at 
zero density the magnetic moment of an electron gas 
is zero by definition. Such approach, however, requires 
the measurements to be performed down to very low 
carrier densities, lower than were accessible earlier. In 
the current paper we extended the thermodynamic mea- 
surements down to densities as low as 0.3 n c , deep in the 
insulating region. This extension enabled us to reveal 
strong unexpected nonlinearity of dm/ dn in weak mag- 
netic fields, which, surprisingly, exists also in the metal 
phase. 

Before presenting the experimental results, we first 
discuss in detail the way we extended the measurements, 
as well as the corrections which should be applied to 
properly interpret the data. Consider a system con- 
sisting of a 2DEG connected through an ohmic contact 
and a battery, which provides the gate voltage V, to the 
gate, as shown in Fig.[T] The free energy of the system 
can be written as follows: 



f = fc + f2D- enV + 



2c ' 



(1) 



where fc and /2D are the free energies of the gate and 
the 2DEG, respectively, and the last term describes the 
electrostatic interaction between the 2DEG and the gate 
electrons. The capacitance cq, defined by Eq. ((TJ), differs 
from the geometric one by a factor of ~ (1 + ^o/dox) -1 ; 
here d ox s=s 190 nm is the oxide thickness for the studied 
Si-MOSFETs, and zq ss 3.5 nm is the average distance 
of the 2D electron layer from the Si/Si02 interface. 

The minimum of / determines the equilibrium den- 
sity of the 2DEG. Differentiating Eq. ([T]) with respect 
to n under the constrain of zero net charge of the sys- 
tem (the variation of the electron density at the gate 
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where /12D and /iq are the chemical potentials of the 
electrons in the 2DEG and the gate respectively. We 
ignored Co dependence on the density; accounting for 
it leads to small corrections for both n, and c defined 
below, of the order of (dco/dn)(n/co). 

A change of an external parameter shifts the equilib- 
rium density n. We assume /ig, the chemical potential 
of a thin Al film, to be independent of B and n. Con- 
sider first the density response on a variation of V, i.e. 
dn/dV = c/e, which defines the gate-to-2DEG capaci- 
tance cllll: 
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Fig. 1. A sketch of the experimental setup, superim- 
posed on the band diagram of a Si MOSFET sample; 
z-direction is perpendicular to the 2DEG plane 

As seen from Eq. (J3J, c is slightly renormalized, by a 
finite 2DEG compressibility. 

Similarly, the density response to magnetic field vari- 
ation is given by: 
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By measuring the recharging current in response to the 
gate voltage modulation one can get the capacitance c, 
and then use it to extract djl/dB. In Eq.[l]we defined jl, 
which contains an additional contribution proportional 
to dco/dB, neglected in [9l [10]. As we shall see below, 
it is indeed small for a 2DEG in Si-MOSFETs. 

Since the thickness of 2DEG is finite, albeit small, an 
in-plane magnetic field couples also to the orbital mo- 
tion, leading to two effects: (i) a diamagnetic shift of the 
spatial quantization levels [HI HH H3] i (h) a change in the 
average distance zq of the electrons from the Si — Si02 
interface due to the asymmetry of the confining poten- 
tial, see Fig.Q] The former effect causes a diamagnetic 
contribution to the chemical potential shift, Sud, which 
should be subtracted, if one is interested in the spin con- 
tribution only: <5/i s = 5[i2D — S^d- The later effect leads 
to a change of the capacitance cq(B), and results in the 
contribution 5[i c , which should be subtracted from 5fl 
in order to get Sfi2D- &U2D = Sjl — Su c with 
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For a non-interacting 2D gas these finite thickness 
(FT) contributions can be found by solving numerically 
the Schrodinger equation in a triangular potential. We 
are interested in the energy of the lowest level |0) of the 
spatial quantization, since it is the only one populated at 
the experimentally relevant densities. We choose Lan- 
dau gauge A = (0,B(z — z ),0), where z = (0\z\0). 
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Fig. 2. Finite thickness contributions to djl/dB at 
B — 1 T; dashed-dotted line is the diamagnetic con- 
tribution dfid/dB, dotted line - dfi c /dB, solid line - 
the sum of the two. Dashed line is dzo/dB (right scale, 
note that zq w 3.5 nm). 

Such a gauge preserves the minimum of the dispersion 
relation eo{k x ,k y ) at k y = 0. Care should be taken 
to average the result over the in-plane k- vectors of the 
populated states to get erj = (eo)fc and zq — (zo)k, since 
components with different k y arc affected differently by 
the field. For density-dependent confining electric field 
E we adopted the expression [TJ 
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where e • no is the depletion region charge, no w 
10 11 cm -2 for our samples. At low, compared to the 
Fermi energy, temperatures the diamagnetic contribu- 
tion dfid/dB = d 2 (ne~o)/dndB. The capacitance con- 
tribution, Eq.[5j can be expressed through zq as: 



d^ c _ 2 d(l/c ) 
dB ~ n€ dB 



Anne 2 dz 
s Sl dB 



(7) 



Both contributions, being the first derivatives of 
even functions, are linear in B at weak magnetic fields; 
they are shown in Fig. [2 For B > the diamag- 
netic contribution is positive, and the capacitance one 
is negative (it is shown in Fig. [5] with negative sign). 
The sum of the FT contributions is positive, there- 
fore disregarding them would lead to underestimation 
of dm s /dn = —d{i s /dB. For our system the resultant 
contribution does not exceed 0.05/is at B = 1 T at low 
densities and drops with n; therefore it is small com- 
pared to the measured dji/dB, at least for low tem- 
peratures, see Fig. |3 The FT contributions should be 
more important for systems with lower effective mass in 
z-direction, such as GaAs based heterostructures, or in 
higher magnetic fields [9 . 

When the gate voltage or the magnetic field is modu- 
lated at a sufficiently small frequency u>, so that the sys- 



tem stays in the thermodynamic equilibrium, the mod- 
ulation leads to a recharging current with the amplitude 

(a) I v = l u J C(u J )SV (b) /b = -^M|§ (5j B, (8) 

e Od 

respectively, where the capacitance C = Sc is propor- 
tional to the sample area S; Eqs. |5] were used in [91 [TU]. 
in order to determine the capacitance (|Ek), and then to 
extract djl/dB (JSJb)- 

However, for low electron densities, recharging of 
the sample is hindered by the large contact and 2DEG 
resistances. When the resistance approaches 1/wC, 
the capacitance becomes complex, and its magnitude 
drops. We show below that, under certain assumptions, 
Eqs. © can still be used for extracting djl/dB using 
the complex capacitance, provided both iy and Ib arc 
measured at the same conditions: frequency, magnetic 
field and temperature. 

Recharging of the capacitor, either by modulation of 
the gate voltage or of the magnetic field, leads to gra- 
dients of the potential and the carrier density over the 
2DEG area, and the electrochemical potential H2D be- 
comes a function of coordinates. In addition, recharging 
current leads to a potential drop across the contact. The 
current density j is governed by the continuity equation: 



dn _^ 

e Tt = ~^ 



V crV 
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where a is, in general, a coordinate-dependent conduc- 
tivity of the 2DEG and the contact region. 

Following [14] we assume that the system stays in 
quasi-thermodynamic equilibrium, and that the char- 
acteristic spatial scale of the electrochemical potential 
variations is large compared to the oxide thickness. 
Then, Eq. © with coordinate-dependent ^2d{t) is sat- 
isfied locally, the assumption called the "local capaci- 
tance approximation" . Under this assumption, accord- 
ing to Eqs. (]3j) and (HJ, both voltage and magnetic field 
modulation generate variation of the density, with the 
only difference being eSV replaced with —dfi/dBSB. 
These variations enter as the source term in Eq. ([9]) , 
whose Green's function, subject to the proper boundary 
conditions, determines the current through the contact. 
This means that Eqs. @ hold under local capacitance 
approximation, with the same C in Eqs. (JS^,) and (JBJa). 

We present results obtained with the most exten- 
sively measured Si-MOSFET sample with the peak mo- 
bility of 3.4m 2 /Vs at 1.7K. This sample is similar to 
that used in Ref. [3]; the metal-insulator transition in 
it occurs at a critical density n c rj 8.5 • 10 10 cm~ 2 . 
The results obtained with many other samples of sim- 
ilar quality were similar. Most of the measurements 
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Fig. 3. djl/dB at n = 5 ■ 10 10 cm -2 , deep in the insu- 
lating regime, full symbols, left scale; the temperatures 
are: ■ - 1.7K • - 3.4K ▲ - 4.6K ▼ - 6.8K. □ - djl/dB 
at n — 2 ■ 10 11 cm -2 and 1.7K, deep in the metallic re- 
gion, right scale. The slope of the dashed line gives the 
derivative of susceptibility —dx/dn, and the slope of 
the solid line - estimated FT contributions, see Fig. [5] 

were done with in-plane magnetic field modulated at fre- 
quency 6.2 Hz with 40 mT amplitude. The results scaled 
linearly with modulation amplitude. Figure [3J shows 
djl/dB as a function of the magnetic field. It is anti- 
symmetric with B at 6.2 Hz, as expected for a derivative 
of an even function. At higher frequencies an asymme- 
try develops, we believe, due to mechanical vibrations; 
it was verified that the antisymmetrized with respect 
to magnetic field djl/dB is frequency-independent up 
to 30 Hz. The upper measurement frequency is limited 
by the EMF in wire loops, as well as by eddy current 
heating. 

From now on let us discuss the results for B > 0. As 
seen in Fig.[3j at the lowest temperature of I.7K, djl/dB 
is strongly nonlinear. For densities deep in the insulat- 
ing regime, e.g. at n = 5 • 10 10 cm~ 2 , as the field in- 
creases, djl/dB sharply drops, then reaches a minimum, 
and saturates, or even increases at larger fields. The 
value of djl/dB = —dm/dn at the minimum is approxi- 
mately —0.7 hb, indicating almost full spin polarization; 
at even lower densities it reaches a value of almost —fiB- 
In the metallic regime, djl/dB changes sign. In Fig. [3] 
we plotted djl/dB for n = 2 • 10 n cm~ 2 , the density at 
which it is maximal. Although the saturation value at 
this density is much smaller, about 0.15 fiB, the satura- 
tion field B* is similar to the one for n = 5 ■ 10 10 cm -2 . 
Note that B* is small, « 0.6 T, well below the polar- 
ization field for a non- interacting system (w 11 T), and 
for the interacting one (rj 2.7 T), both estimated for 
n = 1 • 10 11 cm -2 [SI [7]- Such small B* implies an exis- 
tence of a relevant energy scale, much smaller than the 
bare, and even interaction renormalized Fermi energy. 



The decrease of \dji/dB\ at the lowest temperature 
above B* can be attributed to the finite thickness con- 
tributions: the slope of the solid line in Fig. [3] reason- 
ably agrees with the predicted FT contributions for this 
density. As temperature grows, the nonlinearity smears 
out, and its onset shifts to higher magnetic fields [T2"] . 

To explain qualitatively the results, one needs to in- 
voke, besides itinerant electrons, states, localized either 
within the potential well |15] , or in its close vicinity. 
These states can be recharged at the frequency of the 
measurements, and contribute to d\xir)/dB, but would 
not contribute to Hall conductivity or to Shubnikov-de 
Haas oscillations. The fact that the zero density gate 
voltage Vo (i.e., the so called "threshold voltage"), de- 
termined from extrapolation of the Hall or Shubnikov- 
de Haas measurements, increases with reduction of elec- 
tron mobility 1 indicates the existence of such states. 

On the other hand, the simplistic two-liquid model, 
in which the non-interacting localized states exist in- 
dependently of the interacting 2DEG, suggested in [15], 
cannot explain such a low value of B* , since the charac- 
teristic field, ksT/fis, for polarization of a single spin 
(or non-interacting non-degenerate gas) at 1.7K is 4.2T, 
much bigger than B* . We therefore believe, that any 
explanation of the presented results must involve inter- 
actions. Note, that the underling mechanism should be 
different from the interactions between localized spins, 
discussed in [3J , since such interactions lead to antiferro- 
magnetic coupling, and slower than 1/T divergency of 
the susceptibility, in contrast with our observations (12) . 
We infer that coupling mediated by itinerant electrons 
should play a role. 

Besides localized states, one can think about 
interaction-enhanced FT contributions. Indeed, the 
estimations above ignored interactions. The interaction 
energy, however, is comparable to the energy difference 
between the levels of spatial quantization in the well; 
e.g., for n = 10 11 cm~ 2 both are about 10 meV. It is, 
therefore, not a priori clear that the interactions can 
be ignored in calculations of zo(B). One can envision 
a scenario, in which spin polarization leads to a change 
in zq. Note that Szq w lA would suffice to contribute 
one Bohr magneton to djl/dB. It would be, however, 
quite a coincidence if such a mechanism conspires to 
contribute djl/dB « — [Lb at the saturation. 

In conclusion, we succeeded in measuring the mag- 
netic field dependence of d^D/dB over a density range 
much wider than was accessible earlier, from metallic 
to deep in the insulating region. The dependence of 
dm/dn on B is found to be strongly nonlinear at low 
temperatures, already in magnetic fields much smaller 
than the characteristic fields expected for the 2D elec- 
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tron gas. These results point to a strong coupling be- 
tween the itinerant and localized states in high mobility 
Si-MOSFETs. 
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